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Abstract 

A class of semi-linear hyperbolic systems in 1+1 dimensions was 
investigated several years ago by Ori and Gorbonos [1]. This class, 
to which we shall refer as flux-conserving systems, exhibits a variety 
of interesting mathematical properties. Here we demonstrate how the 
formalism of flux-conserving systems can be applied to the problem of 
black-hole evaporation in 1+1-dimensions. More specifically, we show 
how the semiclassical CGHS [2] field equations may be approximated 
by a certain flux-conserving system. 

1 Introduction 

The investigation of gravitational collapse and black-hole (BE) dynamics 
often leads to systems of partial differential equations (PDEs) in 1+1 di- 
mensions. This happens, for example, in the study of four-dimensional (4D) 
BHs when spherical symmetry is assumed. The usage of spherical symme- 
try for dimensional reduction in BH physics has been a traditional research 
strategy in General Relativity for many decades — since the discovery of the 
Schwarzschild and Reissner- Nordstrom solutions almost a century ago. 

About two decades ago Callan, Giddings, Harvey and Strominger (CGHS) 
[2] introduced a two-dimensional toy model for investigating the semiclassi- 
cal evaporation of black holes, with the hope of addressing the information 
puzzle. For appropriate initial conditions, this model describes the formation 
of a BH as well as its subsequent evaporation. Mathematically, this model 
again reduces to a system of PDEs in 1+1 dimensions. 

Among the family of 4D spherically-symmetric solutions, an interesting 
subclass is the one in which the matter content is made of (at most) two 
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radial fluxes of nuU fluids. The null fluxes are assumed to be pressure-less 
and non-interacting (except gravitationally) , hence each flux is conserved 
separately. If only one such null fluid is present, the system admits an exact 
solution, known as the Vaidya solution |i3j. Obviously, if both fluxes vanish 
the system reduces to the Schwarzschild solution. It is useful to extend this 
subclass by allowing the presence of a spherically-symmetric electric fleld, in 
addition to the null fluids. This gives rise to the mass-inflation model 0], 
the charged Vaidya solution [5j (in the case where one flux vanishes), and, 
obviously, to the Reissner-Nordstrom solution (in the special case where 
both fluxes vanish). All these solutions played a central role in studying 
the non-trivial internal structure of spherically-symmetric charged BHs [6], 
and particularly in investigating the singularity which develops at the inner 
horizon [3 HIS]. 

A few years ago Ori and Gorbonos pj explored a simple class of semi- 
linear hyperbolic systems in 1+1 dimensions, to which I will refer as flux- 
conserving systems. Each such system is characterized by a single function of 
one variable, denoted by h{R) below. This class, which was first introduced 
in Ref. [9j, exhibits several interesting mathematical properties. In partic- 
ular, in any such system one can identify two conserved null fluxes. In the 
special case where one of these fluxes vanish, the solution takes a Vaidya-like 
form and the system then reduces to an ordinary differential equation (ODE), 
which may remarkably simplify the analysis. The flux-conserving systems 
contain, as special cases, the conventional Vaidya [3j and mass-inflation [3] 
solutions (as well as the Schwarzschild and Reissner-Nordstrom solutions), 
which all correspond to the same function h{R) It also includes the 
classical CGHS solution, though the semiclassically-corrected CGHS system 
does not belong to this class. 

The original motivation for introducing these semi-linear systems in Ref. 
[1] was to provide a simpler (and at the same time wider) mathematical 
context for investigating the formation of null weak singularities at the inner 
horizon of BHs. This goal was undertaken by Gorbonos and Wolansky [IQj, 
who demonstrated that for a certain (piecewise-linear) class of functions h{R) 
a generic null weak singularity indeed forms. 

The present paper will focus on a different aspect of these semi-linear 
systems: the very conservation of the individual fluxes, and its possible ap- 
plication to analyzing the evolution of dynamical BHs. As an example, con- 
sider a flux-conserving system which describes a BH bombarded by ingoing 
radiation (with no outflux). In such a system the influx at past null inflnity 
directly determines the rate of change of the BH mass. In the traditional 
framework of conventional 4D spherically-symmetric BHs this is illustrated 
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by the ingoing Vaidya solution, in which the information is fully encoded 
in the mass-function M{v). The same applies to (Vaidya-like solutions for) 
evolving BHs in any type of flux-conserving system. 

More specifically, our main goal in this paper is to demonstrate how the 
mathematical concept of flux-conserving systems may be used to analyze the 
evaporation of two-dimensional semiclassical BHs. 

As was already mentioned above, the classical CGHS system is a flux- 
conserving one [Ij, but this situation changes when the semiclassical effects 
are included. At first glance, an evaporating BH (say, in the Hartle-Hawking 
state) may actually be regarded as an archetype of a system which violates 
flux-conservation: No influx arrives from past null infinity, and yet negative 
influx penetrates into the BH and causes its evaporation. We shall show here, 
however, that by a re-definition of the field variables, the CGHS field equa- 
tions may be transformed into a form which is approximately flux-conserving. 
More specifically, as long as the evaporating BH is still macroscopic, one may 
expand the semiclassical field equations in a small quantity (proportional to 
the number of quantum fields and inversely proportional to the variable R 
defined below). The zero order in this expansion yields the classical system, 
whereas the next order represents the leading semiclassical effects. We shall 
show that up to (and including) this first order, in our new variables the 
transformed field equations constitute a flux-conserving system. This pro- 
vides a powerful analytic tool for analyzing the structure of the semiclassical 
spacetime of a CGHS evaporating BH. 

We start with a very brief review of the formalism of flux-conserving 
systems (Sec. H) and the CGHS model (Sec. HI) — introducing our own 
notation for the latter, and focusing on the properties which are relevant 
to the present analysis. Then in Sec IV we perform the field redefinitions, 
which eliminate the first-order derivative terms from the evolution equa- 
tions. Finally in Sec. V we explore the large- i? approximate system (which 
incorporates the leading-order semiclassical effects), and demonstrate that 
it constitutes a flux-conserving system. We then describe how the CGHS 
evaporating-BH spacetime may be obtained (approximately) from an ingo- 
ing Vaidya-like solution of this system. 

2 Flux conserving systems 

By flux- conserving system we refer to a semi-linear hyperbolic system in 1 + 1 
dimensions, for two unknowns R{u,v) and S{u,v), which takes the form 

Ruv = e^F{R) , S,uv = e^F'{R), (1) 
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for any given function F{R), where a prime denotes d/dR. 

This class of hyperbohc systems displays several unique mathematical 
properties, which were discussed in some detail in Ref. [T]. Here we shall 
briefly mention some of these properties. 

First, the system admits a gauge symmetry: If one transforms the coor- 
dinates as n — )• u{u),v — t- v{v), then R behaves as a scalar (it is unchanged), 
but S transforms as 

(that is, e*^ transforms like a metric component guv) 
Second, we define the two flux functions 

where a comma denotes partial dervative. It then follows from the hyperbolic 
system ([T]) that 

Tuu,v — — 0' (4) 

Namely, T„„ only depends on u, and T^^ anly depends of v. We refer to this 
property as conservation of fluxes. 

Naturally the solutions of this system may be divided into three cate- 
gories: (i) solutions with T^u = Tyy = 0, (ii) solutions in which only one flux 
vanishes; and (iii) solutions in which both fluxes are non-vanishing. We refer 
to case (i) as vacuum-like, and to case (ii) as Vidya-like. Case (iii) is the 
mathematical generalization of the 4D spherically-symmetric systems with 
two null fluids, for example the one studied in Ref. [4]. 

Next we define the "generating function" h{R) by 

F{R) = -h'{R). (5) 

More precisely this is a one-parameter family of functions, due to the arbi- 
trary integration constant. It will often be convenient to express this free 
parameter explicitly in the form 

h{R) = ho{R) - fh, (6) 

where ho{R) denotes a specific integral of ([5]), and rh is an arbitrary constant. 



^The field variables are denoted R, S therein, instead of R, S. Here we use the over-bar 
variables to distinguish them from the specific variables of the CGHS model presented 
below. Similarly, we use T below to distinguish it from the physical stress-energy tensor 
T, and the same for the mass parameter m and the mass function M{v). 
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In the vacuum-like case (i) the space of solutions reduces to a one- 
parameter family of solutions (apart from gauge transformations), character- 
ized by their "mass parameter" fh. It is convenient to express the vacuum- 
like solutions in an "Eddington-like" gauge, in which both R and S depend 
on X = V — u solely. For given fh, R{x) is then determined from the ODE 

^ = ho{R)-rn, (7) 
ax 

and then S{x) by = Hq^R) — fh. The Eddington-like gauge thus exhibits 
the static nature of the vacuum- like solutions. □ 

Assume now that h{R) is monotonically increasing in some range 
Then, at least for some range of the parameter fh, h will vanish at a certain 
^- value denoted ^^nd will become negative at R < Rq. As was discussed 
in Ref. [T], the surface R = Rq functions as an event horizon, and the 
vacuum-like solution then corresponds to a BH. [If h{R) has an additional 
root at R < Rq, it will correspond to an inner horizon.] 

2.1 The Vaidya-like solution 

Assume now that only one flux vanishes. For concreteness we shall consider 
here the ingoing Vaidya-like solution, in which T^^ is nonvanishing. Then, 
as was shown in Ref. [T], each solution is characterized by a mass-function 
M{v). The unknown R{u,v) is determined by the ODE 

R, = ho{R)-M{v). (8) 

In turn the other unknown S{u, v) is given by 

= -Ru- (9) 

The nonvanishing influx satisfles T^^ = M ^ . Note that the solution expressed 
in this from assumes a specific T;-gauge (to which we may refer as Eddington- 
like), but the gauge freedom in u is preserved. 

We shall now derive an expression for S^v which was not included in Ref. 
[1]. Differentiating S from Eq. ([9]) yields S^^ = R^uv/R,u- Substituting now 
Eq. ([1]) for R^uv and Eq. ([9]) again for R^u^ we obtain the simple relation 

5, = -F{R). (10) 



^In cases where the vacuum-hke solution corresponds to a "black hole" (see below), this 
staticity characterizes the external region. Inside the horizon, the vacuum-like solution is 
homogeneous rather than static. The functions R, S then depend on t = v -\- u. 

^This may always be arranged, as long as h{R) is not constant, by transforming R — > 
— _R if necessary (while fixing h). 
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This relation may be useful for constructing approximate solutions for evap- 
orating 1 + 1-dimensional black holes. 

The outgoing Vaidya-like solution — the analogous situation in which the 
nonvanishing flux is rather than T^^ — proceeds in a similar manner, 
though with certain sign changes. One obtains in this case the expressions 
Ru = -[ho{R) - M{u)], = R,, f^u = -M,u, and 5„ = F{R). 



3 The CGHS model 

The CGHS model [2j consists of gravity in 1+1-dimensions coupled to a 
dilaton field (p and to a large number N of identical massless scalar fields 
fi. CGHS provided an action which includes an effective semiclassical term 
derived from the trace anomaly. With appropriate initial conditions, this 
model describes the formation of a BH by a collapsing shell, as well as its 
subsequent evaporation. The field equations are most conveniently expressed 
in double-null coordinates n, v, in which the metric takes the form ds'^ = 
—e^fdudv. The system thus consists of two unknowns, (j)(u,v) and p(u,v). 
The field equations for (p and p are derived in Ref. p]. We find it convenient 
to transform to new variables R = e~'^'^ and S = 2{p — (/)) (as in Ref. 
The evolution equations then take the form 

Rtuv — 6 K p^ny , S,ii^ — K p,ny j R, (If) 

where K = N/12. E The Energy-momentum fiuxes T^u and T^^y are given 
by 

T^uu — R,uS^u R,uu 5 = R,vS^v R,vv (-^■^) 

The system ([TT]) may look quite simple at first glance, but one should 
recall that p = (S — In R)/2 must be substituted in order to close the system, 
which messes the right-hand side. Bringing the equations to their standard 
form, in which R,uv and S,uv are explicitly given in terms of lower-order 



*The above relations for both the ingoing and outgoing Vaidya-like solutions hold 
outside the BH. In the internal region certain sign changes apply, like in the vacuum-like 
case considered above. 

^We have set here the cosmological constant A = 1 by a trivial shift in p, as explained 
in Ref. [11] . Also we set the null solution fi — for the scalar fields, because we are 
interested here in the domain of evaporation where no incoming matter waves are present. 
The system also contains two constraint equations, but they are not needed here because 
any solution of the evolution equations, with appropriate initial conditions, automatically 
satisfies the constraint equations. 
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derivatives, we end up with the system 



_ , i2R-K) K 

K K 

S,uv = (R-K)^ 2i?2 (R-K)- ^'^^^ 

The classical model is obtained by setting K = 0, after which the evo- 
lution equations reduce to R,uv = S,uv = 0, which are easily solved. In 
an Eddington-like gauge the general solution takes the form 

R = e(--) + /„(n) + Mv) , S = v-u, (15) 

where fuiu), fv{v) are arbitrary functions. 

Note that the classical model is a flux-conserving system, with R = 
R, S = S, and ho = R. As such it admits the trivial (vacuum-like) static 
BH solution R = e^'"~'^^ + fh, S = v — u, as well as Vaidya-like solutions. 
This straightforward flux-conserving form breaks down once the semiclassical 
terms are included. 



4 Approximate flux-conserving formulation of CGHS 
dynamics 

4.1 Field redefinition 

The presence of terms quadratic in first-order derivatives makes the field 
equations (|13|14p harder to analyze. We wish to transform the field vari- 
ables so as to get rid of these terms. A straightforward algebra leads to the 
following transformed variables: [12j 

W{R) = ^R{R- K) - Kln{^/R + VR- K) + K (^^ + ln2^ (16) 

and 

Z = S + AZ{R), (17) 

where 

AZ{R) = ^{R-W) - In R. (18) 
K 

We shall only be concerned here about the domain R > K (this includes 
the entire domain-of-dependence portion of the spacetime, both outside 
and inside the BH, because the singularity is located at R = K 
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Since dW/dR = y^l — K/R, the function W{R) is monotonically increasing 
throughout this domain. Therefore, the inverse function R{W) is mathe- 
maticahy weh defined, though it cannot be expressed exphcitly. 

With these new variables the evolution equations take the schematically- 
simpler potential form: 

W^uv = Vw{W) , Z,uv = Vz{W), (19) 

with certain "potentials" Vvi/(^) and VziW). But the simplification does not 
come without a cost: These potentials are explicitly obtained as functions 
of R rather than W . One finds 

Vw = - - -^/L e-^^(^) (20) 

^/R{R-K) 



and 

e-^^(^). (21) 



^ K 



R-K/2 ^ 



Nevertheless, the potential form (jl9p is useful for various analytical consid- 
erations, despite its implicit character. Furthermore, we shall primarily be 
concerned here about the leading-order semiclassical effects, and as it turns 
out, these leading-order effects are fully encoded in the large- (or, equiva- 
lently, large- VF) asymptotic form of the field equations. At the desired order 
the functions VV(l^) and Vz{W) [as well as i2(VF)] all take simple explicit 
forms, as we shortly demonstrate. 



4.2 Large-i? asymptotic behavior 

The semiclassical effects, responsible for the BH evaporation, are encoded in 
the K terms in the field equations (|13ll4p or (|16ti21[) . (Notice that W and 
Z respectively approach R and S as K ^ 0, and the field equations then 
reduce to their classical form, namely Vw = —l,Vz = 0.) To depict the 
leading semiclassical effect we may formally decompose all relevant functions 
in powers of K. Doing so, one observes that this effectively becomes an 
expansion in K/R (see below). Although K itself may be a large number, 
as long as the BH is macroscopic (namely, its mass is ^ K), K/R ^ 1 is 
satisfied throughout the BH exterior (and also in the interior, except in the 
neighborhood of the R = K singularity — a region which will not concern us 
here). Thus, to depict the leading-order semiclassical effect we decompose 
all relevant quantities to first order in K/R (or, equivalently, K/W). For 
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the transformation of variable we obtain 



and 



The inverse function R{W) is given at this order by 



R = W 



K f K\ 

1 + lnT^ + — 

2W \W 



(22) 



(23) 



(24) 



The potentials Vvf(W^)) VziW) can now easily be expanded to first order in 
K/W, 



^ w 



m \w ) 



(25) 



Recall that this leading-order semiclassical approximation is applicable as 
long as KjR ~ KjW <C 1, which for a macroscopic BH holds everywhere in 
the exterior (and also throughout most of the BH interior). 



4.3 The large- i? flux-conserving system 

We now proceed with the above large-i? approximation, and analyze the 
PDE system 

Z,,,, = e^VyXW) (26) 



with 



(27) 



omitting all higher-order terms oc {K/W)"^ . 

Since Vz = dVw/dW, this is a flux-conserving system: It is described by 
Eq. ([1]), with R = W , S = Z , and F = Vw- Its generating function is 



ho{W) = W + ^lnW. 



(28) 



In particular the system (|26l27p admits Vaidya-like solutions. As it turns 
out |13) . at the leading order an evaporating BH may be approximated by 
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an ingoing Vaidya-like solution with a linearly-decreasing mass function. H 
Assuming that the BH was created by the collapse of a shell of mass Mq, 
which propagated along the ingoing null orbit v = vq, the mass function 
takes the form M{v) = Mq — {K/4){v — vq). By shifting the origin of v such 
that vq = —{4:/K)Mq, we obtain 

K 

M{v) = -—V = m^{v). (29) 

Note that dm^/dv = —K/A = —N/AS, which reflects the standard evapora- 
tion rate of a two-dimensional BH [2] . 

The determination of the evaporating-BH spacetime thus reduces within 
this approximation to analyzing the ODE ([8]) which now reads 

W^y = ho{W)-M{v) = {W + ^\nW) + ^v. (30) 



The other unknown Z is then obtained from Eq. (jlOp . which reads 

K 

Zy = l + . (31) 

[Alternatively Z may be obtained from Eq. ([9]).] The dependence upon u 
enters through the initial conditions at v = vq (the collapsing shell). The 
appropriate initial conditions at the shell will be discussed in Ref. [13j. 

Let us denote the solution to the two ODEs (|30|3ip by Wc{u, v) and 
Zc{u,v). Although it constitutes an exact solution to the system (|26l27p . as 
it stands it does not yield an adequately precise description of the CGHS 
BH spacetime. Indeed an approximation was involved in the first place in 
deriving the flux-conserving system (|26|27p from the original CGHS system. 
This approximation is fine, because the terms omitted were of order (K/W)'^, 
whereas our attempt here is to construct an approximate solution at order 
K/W. As it turns out [13), however, the solution {Wc,Zc) fails to satisfy 
the proper initial conditions at the shell {v = vq), and the discrepancy is 
at order K/W, which cannot be neglected here. Stated in other words, 
the flux-conserving solution which properly describes the evaporating BH is 
not a pure ingoing Vaidya-like solution, but one which also contains a weak 
outflux — namely, a small correction 5W{u) (x K to W = Wc- Nevertheless, 
owing to the smallness of 6W (and to the trivial form of the field equations 
at the limit if — )• 0), one can simply superpose the correction 6W{u) on top 



^However, a weak outflux term oc (K/W) has to be added in order to match the initial 
conditions, as discussed below. 
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of the Vaidya-like solution Wc (with Zc unchanged) |13| . This procedure 
yields a sufficiently accurate solution which satisfies the appropriate initial 
conditions. 

In a forthcoming paper jl3j we shall implement this formalism, and con- 
struct an approximate solution R{u,v), S{u,v) to the CGHS field equations 
which satisfies the correct initial conditions (at the shell as well as at past 
null infinity), and is accurate at first order in K/W. 
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